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摘要 : 本 文 提出 了 一 个 数学 启发 算法 求解 单 源 设施 区 位 问题 (SSCFLP)。 该 算法 从 
个 初始 解 开 始 ， 迭 代 地 对 当前 解 进 行 超 大 邻 域 搜 索 改 进 ， 直 到 若干 次 尝试 不 能 该 改进 当 
前 解 为 止 。 算 法 中 ， 初 始 解 使 用 拉 格 朗 日 松弛 启发 算法 ， 超 大 邻 域 搜索 采用 子 问题 数学 
模型 精确 求解 。 算 法 设计 的 要 点 在 于 超大 邻 域 的 选择 , 既 不 能 太 大 造成 子 问 题 求解 困难 ， 
也 不 能 太 小 造成 当前 解难 以 改进 。 使 用 5 组 272 个 基准 测试 案例 进行 算法 测试 ， 共 发 现 
191 个 案例 的 最 优 解 ， 更 新 36 个 案例 的 已 知 最 好 解 ， 表 明 数 学 启发 算法 性 能 优异 。 与 近 
年 的 代表 性 算法 (如 制 平面 、 核 搜索 、 超 图 多 交换 启发 、 亡 道 方法 ) 进 行 比较 ， 在 大 规模 
案例 上 ， 本 文 算法 无 论 求解 质量 还 是 计算 时 间 均 具有 显著 的 优势 。 

关键 词 : 单 源 设施 区 位 问题 数学 启发 算法 ; 超大 规模 邻 域 搜索 ; 子 问题 
引言 

设施 区 位 问题 (FLP) 是 一 大 类 寻找 最 优 设施 位 置 的 问题 。 该 问题 广泛 应 用 于 学 校 、 医 
疗 、 公 安 、 养 老 、 应 急 、 物 流 等 公共 或 商业 设施 选 址 决策 。 根据 应 用 场景 的 不 同 ，FLP 类 
型 众多 : 连续 空间 区 位 或 离散 区 位 、 有 无 设施 容量 限制 、 需 求 是 否 允 许 拆 分 指派 、 是 否 
考虑 设施 成 本 、 是 否 指定 设施 数量 等 ， 规 划 目 标 划 分 为 效率 目标 、 公 平 目标 等 。 根 据 设 
施 服务 与 需求 是 否 存在 不 确定 性 ， 又 可 将 FLP 划分 为 确定 性 FLP 和 随机 性 FLP。 其 中 ， 
单 源 CFLP 问题 (SSCFLP) 是 需求 不 允许 拆 分 指派 的 、 有 容量 约束 的 、 顾 及 设施 成 本 的 离散 
区 位 问题 ， 求 解难 度 较 高 ， 也 是 近年 关注 度 最 高 的 区 位 问题 之 一 。 

令 集合 /为 潜在 设施 区 位 ， 和 集合 ./ 为 客户 区 位 ， 设 施 i CED 有 最 大 服务 容量 w 和 固定 
建设 成 本 为 fp， 客 户 j (I E 力 有 服务 需求 4;， 被 设施 i 服务 的 费用 为 cy。 定 义 布尔 型 决策 变 
量 yj 为 是 否 在 区 位 ;建设 施 ， 布 尔 型 决策 变量 wj 为 是 否 指派 设施 i 服务 客户 j，SSCFLP 数 学 


模型 如 下 : 
Minimize  Yier fi + Lier Djey CijXij (1) 
Subject to: ier xij = 1, Vj EJ (2) 
È jej djXij < SiYi, Vi EI (3) 
yi = {0,1}, viel (4) 
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xij = {0,1} Viel jes (5) 
模型 中 , 目标 函数 (1) 最 小 化 设施 开设 成 本 和 获取 服务 成 本 。 约束 (2) 确 保 每 个 客户 由 唯一 
的 设施 提供 服务 ; 约束 (3) 要 求 指派 给 每 个 设施 的 客户 需求 量 不 超过 其 最 大 服务 容量 。 的 
束 (4) 和 (5) 定 义 两 类 布尔 型 决策 变量 。 

求解 SSCFLP 问题 的 方法 划分 为 两 大 类 : 精确 算法 和 启发 式 算 法 (Basu et al. 2015, 
Ulukan & Demircioğlu 2015)。 精 确 算法 包括 分 枝 定 界 (Neebe & Rao 1983, Holmberg et al. 
1999)、 列 生成 (Dfaz & Fernandez 2002)、Cut-and-Solve (Yang et al. 2012, Gadegaard et al. 
2018) 等 。 因 SSCFLP 问题 具有 极 高 的 计算 复杂 度 , 精确 算法 难以 求解 规模 较 大 的 案例 。 

与 精确 算法 相 比 ， 求 解 SSCFLP 的 启发 式 算法 相当 多 样 。 基 于 拉 格 朗 日 松弛 的 启发 
式 算 法 (LH) 应 用 最 广泛 (Barcelo & Casanova 1984; Klincewicz & Luss 1986; Beasley 1993; 
Sridharan 1993; Agar & Salhi 1998; Hindi & Pienkosz 1999; Cortinhal & Captivo 2003). 1% 
方法 将 指派 约束 或 容量 约束 整合 到 目标 函数 中 ， 然 后 迭代 执行 以 下 步骤 : 中 对 偶 问题 求 
解 ， 更 新 原 问题 下 界 ; @ 利 用 对 侦 问 题解 构造 原 问 题 可 行 解 ， 更 新 原 问 题 上 界 ; OM 
度 方 法 更 新 乘 子 。 基 于 提供 的 问题 下 界 ，LH 算法 常常 和 其 他 算法 整合 使 用 。 例 如 ， 
Holmberg et al. (1999) 将 LH 算法 与 分 枝 定 界 精确 算法 框架 整合 ，Chen & Ting (2008) 将 
LH 算法 与 蚁 群 系统 混合 。 其 他 启发 式 算法 包括 : 禁忌 搜索 (Filho & Galvaol998; Delmaire 
et al. 1999; Cortinhal & Captivo 2003)、 超 大 规模 邻 域 搜索 (VLNS) (Ahuja et al. 2004; Tran 
et al. 2017). EIER (Contreras & Diaz, 2008 )、 核 搜索 (Kernel Search)(Guastaroba & 
Speranza 2014)、“ 廊 道 " 方 法 (Corridor Method)(Caserta & Voß, 2020) 等 。 

近 十 余年 来 , 求解 SSCFLP 的 代表 性 算法 如 下 : 核 搜 索 (Guastaroba & Speranza 2014), 
割 平 面 (Yang et al. 2012; Gadegaard et al. 2018)、 超 大 邻 域 搜 索 (Tran et al. 2017) 和 数学 启 
发 (Caserta & Voß, 2020)。 这 些 算 法 在 求解 质量 方面 具有 竞争 力 ， 也 测试 了 多 组 案例 ， 验 
证 了 算法 的 通用 性 。 因 这 些 算法 计算 时 间 偏 长 ，Oliveira et al.(2020) 设计 了 松弛 自 适 应 
记忆 规划 (RAMP) 方 法 ， 大 幅 降 低 计 算 时 间 ， 获 得 较为 满意 的 结果 。 

本 文 提 出 一 个 数学 启发 式 算法 求解 SSCFLP 案例 。 五 组 272 个 基准 测试 案例 的 实验 
表明 ， 数 学 启发 式 算法 在 求解 大 规模 SSCFLP 案例 中 具有 竞争 优势 : 计算 结果 优 于 目前 
公认 最 好 的 几 个 算法 ， 且 求解 效率 显著 提升 。 

1 数学 启发 算法 

邻 域 搜索 用 于 渐进 改进 当前 解 , 是 众多 优化 问题 启发 算法 的 关键 步骤 .区 位 问题 中 ， 
设施 的 变动 必然 会 引起 需求 指派 变化 ， 使 得 邻 域 搜索 较为 复杂 。 因 此 ， 区 位 问题 算法 设 
计 的 一 个 核心 问题 是 : 如 何 有 效 地 移动 设施 并 更 新 指派 ， 寻 找到 可 行 的 、 更 节约 的 邻 域 
解 。 考 虑 到 SSCFLP 存在 设施 容量 约束 ， 涉 及 较 多 设施 的 超大 邻 域 搜 索 才 有 较 高 可 能 性 
改进 当前 解 。 然 而 ，SSCFLP 的 大 邻 域 结构 定义 复杂 ， 邻 域 搜索 的 计算 复杂 度 也 很 高 
(Ahuja et al. 2004; Tran et al. 2017)， 均 存在 指数 级 数量 的 客户 移动 移动 和 设施 移动 。 

本 文 拟 采用 数学 模型 求解 SSCFLP 的 超大 邻 域 。 算 法 原理 如 下 : OEH LH 方法 
(Holmberg et al. 1999) 获 得 问题 初始 解 ; 包 从 当前 解 中 随机 选择 一 个 超大 邻 域 ， 获 得 邻 域 


内 需求 点 、 当 前 设施 和 候选 设施 集合 ; OMET eA IF, E 
指派 ; 外 重复 执行 步骤 @ 和 人 @， 直 到 若干 次 尝试 不 能 更 新 当前 解 。 步 


新 设施 集合 和 需求 
PROMI — TAB 


域 ， 步 又 G9) 用 于 发 现 邻 域内 的 最 优 解 。 基 于 这 一 思路 ， 数 学 启发 算法 示意 如 下 : 
参数 :连续 未 更 新 最 好 解 循环 数 (mloops) 


1. spest = LagrangianHeuristic(psize); 


2. notimpr=0; 


3. While notimpr < mloops: 


4. 
5, 
6. 
7. 
8. 


Select a large neighborhood /* and J* from solution s; 
s*= SolveSubproblem(/*, J*); 
s’=UpdateSolution(s,s*) 
If f(5’)<f(Sbest): Spest=s’, notImpr=0; 
else: notimpr+=1; 


9. Output Sbest- 


Ft 


y 法 中 ， 步 又 (1) 使 用 拉 格 朗 日 松弛 启发 方法 (Holmberg et al. 1999) 获 得 初始 解 。 


以 上 


F 


模型 步骤 (4) 从 当前 解 中 随机 选择 一 个 超大 邻 域 , 包括 设施 子 集 /* 和 需求 子 集 J*。 步 又 (5) 


zm 


使 


用 设施 子 集 和 需求 子 集 构造 子 问题 模型 ， 使 用 MIP 优化 器 求解 ， 获 得 局 部 解 s*。 步 又 


(6) 使 用 局 部 解 更 新 当前 解 ， 获 得 新 解 5。 步 骤 (7) 尝 试 更 新 最 好 解 。 当 若干 次 循环 未 更 新 


最 好 解 ， 


从 当前 解 中 选择 一 个 超大 邻 域 (步骤 4) 是 算法 设计 的 要 点 。 邻 域 选择 步骤 如 下 : Ob 
机 选择 一 个 需求 点 ， 以 该 点 为 中 心 ， 从 当前 设施 集合 中 自 近 而 远 地 选 择 K 个 ; QIK 
K 个 设施 的 需求 点 构成 需求 点 子 集 J*; @@) 子 集 J* 中 每 个 需求 点 最 近 的 设施 以 及 已 选择 的 
K 个 设施 构成 子 集 !*。 以 上 步 又 所 选择 令 域 的 大 小 受 K 值 影响 ，K 值 取 值 范围 为 [Kmi， 
Kmax]， 其 中 ，Kmin= min(L/2,7)，Kmox=min(L,10),， 上 为 当前 设施 数量 。 当 L<7 时 , K 处 于 区 


os 
FF 


法 终止。 


间 [L/2,L]; 当 8<L<10 IN, K 处 于 区 则 [L/2,4]; “4 11<L<13 时 ，K 处 于 区 间 [L/2,10]; “4 L214 
HT, K 处 于 区 间 [7,10]。 若 邻 域 太 小 ， 设 施 区 位 往往 难以 更 新 ， 若 邻 域 过 大 ， 模 型 计算 效 


率 低 。 


以 上 邻 域 选择 K 的 取 值 可 以 保证 邻 域 足 够 大 , 但 可 能 会 使 模型 计算 效率 低下 。 为 此 ， 
当 集 合 /* 的 规模 偏 大 时 ， 需 要 精简 一 部 分 候选 设施 ， 再 进行 子 问题 建 模 和 求解 。 精 简 候 
选 设 施 采用 以 下 步骤 : 第 一 ， 若 当前 解 目标 值 小 于 某 个 设施 对 应 的 排除 姜 值 ， 则 这 个 设 
施 不 在 最 优 解 中 。 设 施 排除 方法 来 自 文献 (Holmberg et al. 1999) 提 供 的 LH 算法 ， 即 针对 
拉 格 朗 日 乘 子 ， 计 算 对 偶 问题 中 每 个 设施 子 问 题目 标 值 和 对 偶 问 题目 标 值 ， 若 某 个 设施 
子 问题 目标 值 与 对 偶 问 题目 标 值 之 和 大 于 原 问 题 上 界 , 则 该 设施 不 在 最 优 解 中 。 本 文 LH 
算法 记录 并 更 新 每 个 设施 对 应 的 被 排除 目标 值 闵 值 ， 随 着 解 的 不 断 改 进 ， 可 能 有 越 来 越 
多 的 设施 被 排除 在 最 优 解 之 外 。 第 二 ， 排 除 部 分 设施 后 ， 如 果 设 施 仍然 比较 多 ， 先 选择 
K 个 已 使 用 设施 ， 再 从 其 余 设 施 中 随机 选择 K 个 设施 。 


2 算法 测试 结果 


2.1 WERKA 


本 文 收集 5 个 SSCFLP 基准 案例 集 ， 包 含 272 个 实例 。 案 例 集 名 称 、 来 源 、 分 组 及 
规模 汇总 于 表 1。 案 例 集 OR-Library、Holmberg ~ Yang 和 Tebdl 数据 见 网 站 https://or- 
brescia. unibs.it/instances/instances_sscflp; 案例 集 TB4 Jl https://github.com /SuneGadegaard 
/SSCFLPsolver。 表 中 , 每 个 案例 集 划分 为 若干 分 组 , 分 组 案例 名 称 列 在 分 组 后 面 括号 内 ; 


|7| 和 || 分 表 表 示 候 选 设施 数量 和 需求 点 数量 。 这 些 案例 集 在 供需 规模 、 分 布 、 比 例 


等 方面 均 具有 代表 性 ， 广 泛 应 用 于 区 位 问题 算法 测试 。 
表 1 基准 测试 案例 集 


Mm 


案例 集 分 组 (案例 名 称 ) 案例 数量 I J 
OR-Lib ORI (cap61-cap74) 8 16 50 
(Ahuja et al. 2004) OR2 (cap91-cap104) 8 25 50 
OR3 (cap121-cap134) 8 50 50 
OR4 (capax, capbx, capex) 12 100 1000 
Holmberg H1 (pl-p12) 12 10 50 
(Holmberg et al. 1999) H2 (p13-p24) 12 20 50 
H3 (p25-p40) 16 30 150 
H4 (p41-p55) 15 10-30 70-100 
H5 (p56-p71) 16 30 200 
Yang Y1 (30_200_x) 5 30 200 
(Yang et al. 2012) Y2 (60_200_x) 5 60 200 
Y3 (60 300 x) 5 60 300 
Y4 (80_400_x) 5 80 400 
TB4 G1 (50_100 x x) 15 50 100 
(Gadegaard et al. 2018) G2 (50_200 x x) 15 50 200 
G3 (60 300 x x) 15 60 300 
Tebd1l T1 (1300 x) 20 300 300 
(Avella and Boccia 2009) T2 (i3001500_x) 20 300 1500 
T3 (1500 x) 20 500 500 
T4 (1700_ x) 20 700 700 
T5 (i1000 x) 20 1000 1000 


2.2 Etr 


本 文 数学 启发 算法 使 用 Python 程序 设计 语言 编程 实现 。 算 法 步骤 如 下 : 读数 据 , 使 
用 LH 方法 获得 一 组 初始 解 ， 友 代 执 行 邻 域 选 择 和 子 问 题 求解 ， 直 到 满足 终止 条 件 。 子 
问题 求解 环节 , 使 用 PuLP(https://github.com/coin-or/pulp) 线 性 规划 建 模 工具 完成 SSCFLP 
子 问 题 模型 建构 ， 再 调用 IBM CPLEX 12.6 优化 器 完成 模型 求解 。 为 提升 计算 速度 ， 算 法 


在 PyPy7 (https://www.pypy.org) 环 境 
Core 17-6700 CPU 3.40-GHz 和 8GB 内 存 ，Windows 10 操作 系统 。 


' 运 行 。 实 验 计 算 环境 为 : HP 桌面 计算 机 , 配置 intel 


案例 测试 过 程 如 下 。 首 先 ， 设置 参数 mloops， 即 连续 未 更 新 最 好 解 循环 数 。 每 组 案 


例 使 用 参数 依据 经 验 设 置 ， 案 例 OR1、OR2、OR3、H1 和 H2 设置 为 10， 案 例 OR4、T2 
设置 为 50， 案 例 H3、H4 和 H5 设置 为 20， 其 原理 设置 为 100。 一般 来 说 ， 参 数 psize 数 
值 越 大 ， 计 算 结果 越 稳定 ， 所 需 计算 时 间 也 越 长 。 其 次 ， 针 对 每 个 案例 ， 使 用 数学 启发 


算法 计算 5 次 ,统计 5 个 目标 值 的 最 小 值 、 平 均值 、 最 大 值 和 标准 差 ， 以 及 平均 计算 时 


间 。 平 均值 与 案例 目标 值 下 界 之 差 反 映 算 法 求解 质量 ， 标 准 差 用 于 评估 算法 的 稳健 性 ， 
计算 时 间 用 于 评价 算法 的 计算 效率 。 


2.3 MRAR 
使 用 本 文 数学 启发 算 


界 的 差异 值 Gap, Gap 值 越 小 说 明 


法 对 5 个 案例 集中 272 个 案例 进行 求解 ， 结 果 统 计 如 表 2。 


Fe 


个 案例 计算 5 次 ， 统 计 目 标 值 平均 值 和 计算 时 间 平 均值 。 再 计算 目标 值 与 案例 目标 值 下 


F- 


miy 


= 


法 求解 质量 越 高 。 当 案例 最 优 解 已 知 时 ，Gap=( 目 标 
直 平 均值 -最 优 解 目 标 值 )/ 最 优 解 目 标 值 *100%; 当 案 例 最 优 解 未 知 时 ，Gap=( 目 标 值 平均 
-目标 值 下 界 )/ 目 标 值 下 界 *100%。 表 中 ， 


pin 


“Exact” 栏 列 出 精确 求解 获得 最 优 解 的 数量 


“#Opt” 以 及 案例 平均 求解 时 间 , 来 自 CPLEX 分 支 切 割 算 法 或 改进 割 平 面 算法 计算 结果 。 


Ja 


“MH” 栏 列 出 本 文 数 学 启发 


为 核 搜索 算法 (Guastaroba & Speranza 2014) 


F- 


法 求解 结果 的 Gap 平均 值 和 计算 时 间 平 均值 ，“KSs” 栏 


结果 统计 ; “HMEH” 栏 为 基于 超 图 的 多 交 


换 启 发 算法 (Tran et al. 2017 ) 求 解 结果 统 计 ;“CM” 栏 为 “ 廊 道 ”方法 (Caserta & Voß 2020) 
结果 统计 。 注 意 ， 表 中 案例 求解 质量 可 以 直接 对 比 ， 但 计算 时 间 不 宜 直接 对 比 。 
表 2 案例 计算 结果 统计 


zo 淹 逢 一 


Gap/% Time/s Gap/% Time/s Gap/% Time/s 


1.22 
2.23 
1.56 
149.06 
1.13 
2.13 
8.64 
3.88 
15.59 
76.31 
27.78 
120.58 
232.79 
86.55 
47.83 
61.59 
60.88 
87.12 
158.69 
341.18 


Exact 
#opt Time/s Gap/% Time/s 
OR-Lib ORI 8/8 0.05 0.00 
OR2 8/8 0.06 0.00 
OR3 8/8 0.08 0.00 
OR4 12/12 112.44 0.01 
Holm. H1 12/12 0.20 0.00 
H2 12/12 0.34 0.02 
H3 16/16 2.61 0.00 
H4 15/15 0.67 0.00 
H5 16/16 5.29 0.00 
Yang Y1 5/5 51.00 0.00 
Y2 5/5 1261.82 0.01 
Y3 5/5 65.63 0.04 
Y4 5/5 228.01 0.09 
TB4 Gl 15/15 676 0.03 
G2 13/15 4036 0.07 
G3 10/15 14617 0.07 
Tebdl T1 13/20 3722 0.15 
T2 20/20 47 0.00 
T3 8/20 2017 0.27 
T4 4/20 7744 0.47 
T5 0/20 8275 0.57 


345.58 


0.00 0.29 - - = - 
0.00 0.39 - - E - 
0.00 0.62 - - 7 - 
0.00 34.67 0.04 42.67 9.00 43 
0.00 0.32 3 = 
0.00 0.42 
0.00 0.38 = = 
0.00 2.43 0.00 4.08 z 7 
0.00 0.54 0.00 1.08 = 3 
0.00 2.32 0.00 15.53 7 - 
0.00 411.28 - - - - 
0.00 1640.42 - - 7 - 
0.00 597.06 - - - - 
0.00 1409.11 - - - - 


0.56 2206.96 0.54 428.03 0.23 807 
0.00 334.71 0.01 1159.33 0.00 29 
0.66 4190.28 0.52 2982.72 0.36 1024 
0.90 5244.69 0.82 4992.14 0.78 912 
1.07 6533.15 1.10 8582.74 1.11 932 


表 2 中 ， 案 例 OR1, OR2 和 OR3 由 本 文 作者 使 用 CPLEX 12.6 求解 ， 获 得 全 部 最 优 
fief; 案例 OR4、T1、T2、T3、T4 和 T5 是 使 用 CPLEX12.2 的 计算 结果 (Guastaroba et al. 


2014); 案例 集 Holmberg. Yang 和 TB4 it 


算 结 果 来 自 改 进 割 平面 算法 (Gadegaard et al. 


2018)。 据 文献 Gadegaard et al. (2018)， 改 进 割 平面 算法 求解 TB4 案例 ， 获 得 41 个 案例 


的 最 优 解 。 作 者 发 现 ， 其 


1 3 个 案例 的 最 优 解 有 错误 ， 因 此 表 1 中 案例 集 TB4 共有 38 


个 最 优 解 。 案 例 集 Tebd1 的 最 优 解 目 标 值 、 最 优 解 案 例 数量 和 计算 时 间 引 自 文献 Caserta 


& VoB (2020)。 表 中 可 见 ， 精 确 算法 能 够 高 效 地 求解 规模 较 小 的 案例 ， 如 大 部 分 OR- 
Library, Holmberg 案例 。 然 而 ， 案 例 集 Yang. TB4 和 Tbed1 中 多 数 案例 计算 时 间 很 长 。 


据 文 献 Yang et al.(2014)， 使 用 CPLEX 求解 Yang 案例 ， 每 个 案例 计算 时 间 限 制 为 50000 


秒 ，20 个 案例 中 14 个 案例 找 不 到 最 优 解 。 据 文献 Caserta & Voß (2020)， 使 用 CPLEX 求 
综 上 ， 精 确 算法 求解 中 大 规模 SSCFLP R 


Tbed1 案例 ，100 个 案例 有 45 个 获得 最 优 解 。 
例 存在 困难 ， 需 要 设计 启发 式 算法 求解 这 一 问 


从 表 2 种 可 以 发 现 : (1) 对 于 较 小 规模 案例 (OR1、OR2 和 OR3)， 启 发 式 算法 中 ， 核 搜 


题 。 


索 算 法 表现 最 好 , 但 CPLEX 分 支 切 割 算 法 求解 更 高 效 。(2) 随 案例 规模 增 大 ,对 于 Holmberg 


和 OR4 案例 , 核 搜索 算法 表现 最 好 , 在 求解 质量 与 计算 


时 间 方 面 与 割 平面 算法 差异 不 明 


显 。(3) 对 于 Yang 案例 ，CPLEX 精确 求解 困难 ， 割 平面 方法 大 幅 降 低 了 求解 时 间 ， 且 获得 
最 优 解 ， 核 搜索 算法 能 够 获得 最 优 解 ， 但 计算 时 间 较 长 ， 数 学 启发 能 够 在 短 时 间 内 获得 


高 质量 解 。(4) 对 于 TB4 案例 ，CPLEX 精确 求解 困难 ， 审 


| 平面 方法 大 幅 降低 了 求解 时 间 ， 


获得 大 多 数 案例 最 优 解 ， 而 数学 启发 能 够 在 短 时 间 内 获得 高 质量 解 。(5) 对 于 更 大 规模 案 
例 集 Tbed1， 割 平面 算法 的 可 行 性 尚未 验证 ，CPLEX 能 够 高 效 求解 部 分 案例 ， 如 T2 案例 


算法 核 搜索 和 大 邻 域 搜索 算法 能 较 好 地 求解 Tbed1 
性 能 更 优 ， 特 别 是 对 于 较 大 规模 案例 T4 和 T5， 求 解 质量 明显 提升 ， 且 计算 时 间 大 幅度 


和 T1 中 少数 案例 ， 但 对 于 其 他 案例 ， 精 确 求 解 效率 偏 低 ， 或 者 难以 精确 求解 ， 局 发 式 


! 的 大 规模 案例 , 本 文 数 学 启发 算法 


下 降 。 综 上 分 析 , 数学 启发 算法 在 求解 大 规模 案例 (如 T3、T4、T5) 或 者 较 困难 案例 (如 TB4) 


方面 具有 明显 的 优势 。 


使 用 数学 启发 算法 进行 案例 测试 过 程 中 , 作者 发 现 现 有 文献 中 存在 某 些 错误 。 首先， 
TB4 案例 集中 ，3 个 案例 (50-200-2-4、50-200-2-5 和 60-200-2-5) 的 目标 值 小 于 文献 
Gadegaard et al. (2018) 提 供 的 最 优 目 标 值 。 使 用 CPLEX 重新 计算 案例 模型 ， 结 果 表 明 该 
文献 提供 的 案例 最 优 解 错 误 。 第 二 ，Tbed1 案例 集中 ， 部 分 案例 目标 值 小 于 文献 


Guastaroba et al. (2014) 提 供 的 最 优 目 标 值 ， 应 采用 文献 Caserta & VoB (2020) 提 供 了 新 的 
CPLEX 计算 结果 。 第 三 ， 文 献 Tran etal. (2017) 更 新 的 最 好 解 中 ， 个 别 案例 目标 值 显著 小 
于 案例 最 优 目 标 值 ， 如 Tbed1 案例 i300_7、i300_10 和 i500_13。 另 外 ， 作 者 使 用 CPLEX 
成 功 求解 个 别 案例 ， 更 新 了 最 优 解 ， 如 Tbed1 案例 i300_6、i300 7、i500_6、i500_7、 
i500 9、i500_10、i500_13、i500_15、i700_14 和 i700_20。 表 3 为 部 分 案例 目标 值 上 界 、 


下 界 和 最 优 值 的 更 新 或 更 正 ， 标 注 星 号 的 为 最 


优 值 。 


表 4 汇 总 常见 算法 发 现 最 优 解 的 数量 。 表 中 ， 案 例 已 知 最 优 解数 量 “#Opt” 来 自 文 


献 和 作者 计算 结果 ， 己 统计 表 3 所 列 案例 。 己 知 最 好 解 引 自 文献 Yang et al. (2012)、 
Guastaroba et al.(2014)、Tran et al. (2017) 和 Gadegaard et al. (2018)。 统 计时 ， 剔 除了 文献 
Tran et al. (2017) 和 Gadegaard et al. (2018) 中 有 错误 的 案 


HMEH、CM 分 别 表 示 数 学 启发 算法 、 割 平面 算 


方法 。 对 于 20 个 Yang 案例 ，CPLEX 求解 困难 ， 
et al. 2018) 能 够 发 现 全 部 案例 的 最 优 解 ， 核 搜索 启发 算 


割 平 面 算 


网。 表 中 ，MH、CS、CS2、KS、 
法 、 改 进 制 平面 算法 、 核 搜索 算法 和 亡 道 
法 (Yang et al. 2012, Gadegaard 


法 和 本 文 数学 


启发 算法 均 发 现 了 


18 个 最 优 解 。 对 于 45 个 TB4 案例 ， 改 进 割 平面 算法 (Gadegaard et al. 2018) 获得 38 个 案 


例 的 最 优 解 (已 剔除 3 个 有 错误 的 案例 )， 本 文 数 学 启发 算法 发 现 33 个 案例 的 最 优 解 、 更 
新 4 个 案例 的 最 好 解 ， 计 算 时 间 大 幅 降低 。 对 于 100 个 大 规模 案例 Tbed1，CPLEX 计算 
结果 共 发 现 55 个 案例 的 最 优 解 ,但 计算 时 间 偏 长 ; 核 搜 索 启 发 算法 获得 40 个 案例 最 优 
解 ，HMEH 算法 发 现 28 个 案例 最 优 解 ， 亡 道 方法 发 现 44 个 最 优 解 ， 本 文 数学 启发 算法 
获得 41 案例 最 优 解 ， 计 算 时 间 大 幅 降低 。 

表 3 部 分 案例 CPLEX 计算 结果 


入 


iu 


案例 集 案例 目标 下 界 ”目标 上 界 备注 
TB4 50-200-2-4 25951.63 25955 IE(Gadegaard et al. 2018) 
TB4 50-200-2-5 25326.5 25329 ”更 正 (Gadegaard et al. 2018) 
TB4 50-300-2-5 37138.1 37142 “更正 (Gadegaard et al. 2018) 
Tbed1 i300 6 11326.43* 11326.43 ”更 新 (Caserta & Voß 2020) 
Tbedl i300 7 11470.31* 11470.31 更 新 (Caserta & Voß 2020) 
Tbed1 i500_6 15853.35* 15853.35 ”更 新 (Caserta & Voß 2020) 
Tbed1 i500_7 16205.15* 16205.15 ”更 新 (Caserta & Voß 2020) 
Tbed1 i500_9 16399.40* ”16399.4 ”更 新 (Caserta & Voß 2020) 
Tbed1 i500_10 15886.54* 15886.54 ”更 新 (Caserta & Voß 2020) 
Tbed1 i500_13 13715.96* 13715.96 更 新 (Caserta & Voß 2020) 
Tbedl i500_15 13947.12* 13947.12 ”更 新 (Caserta & Voß 2020) 
Tbed1 i700_14 17383.87* 17383.87 ”更 新 (Caserta & Voß 2020) 
Tbed1 i700_20 15492.02* 15492.02 ”更 新 (Caserta & Voß 2020) 
表 4 常见 算法 发 现 最 优 解 的 案例 数量 


案例 集 分 组 # 案 例 #Opt CPLEX MH CS CS2 KS HMEH CM 


OR-Lib OR4 12 12 12 11 - - 12 - 12 
Yang Yl 5 5 3 5 5 5 5 - - 
Y2 5 5 1 4 5 5 4 - - 
Y3 5 5 1 5 5 5 4 - - 
Y4 5 5 1 4 5 5 5 - - 
TB4 G1 15 15 - 12 - 15 - - 
G2 15 13 - 7 - 13 - - - 
G3 15 10 - 7 - 10 - - - 
Tebd1 T1 20 15 15 12 - - 12 10 12 
T2 20 20 20 17 - - 20 14 20 
T3 20 14 14 7 - - 6 4 8 
T4 20 6 6 5 - - 2 0 4 
T5 20 0 0 0 - - 0 0 0 


对 于 53 个 最 优 解 未 知 的 TB4 和 Tbedl 案例 ， 本 文 更 新 了 其 中 36 个 案例 的 已 知 最 
好 解 ， 见 表 5。 表 中 包括 4 个 TB4 案例 和 32 个 Tbedl 案例 。TB4 案例 目标 值 下 界 和 已 
知 最 好 解 目 标 值 由 改进 割 平面 算法 (Gadegaard et al. 2018)) 提 供 ; Tbedl 案例 目标 值 下 界 
来 自 CPLEX 求解 ， 已 知 最 好 解 目 标 值 来 自 CPLEX、 核 搜索 、HMEH 或 廊 道 方法 。 新 最 
好 解 是 本 文 算 法 求解 结果 。Tbedl 案例 的 最 好 解 处 于 持续 更 新 中 ，HMEH 算法 、 廊 道 算 
法 更 新 了 部 分 案例 最 好 解 ， 本 文 算法 进一步 更 新 了 其 中 32 个 案例 的 最 好 解 。 与 原 有 最 


好 解 相 比 , 这 些 案例 目标 值 平 均 改 进 了 0.48%. 其 中 , 案例 60-300-2-4. i700_1. i700_2. 
i700 4, i1000 1、il000 2、il000 3、il000 4, i1000 5 和 il000 10 的 改进 幅度 明显 ， 
处 于 0.75%~3.55% 区 间 。 


表 5 更 新 已 知 最 好 解 目 标 值 


案例 集 ”案例 下 界 BAREM ”新 最 好 解 

TB4 60-300-2-1 34858.5 34861 34860 
TB4 60-300-2-2 36543.5 36742 36551 
TB4 60-300-2-3 34876.2 34884 34878 
TB4 60-300-2-4 34817.6 36057 34821 
Tbed1 1300_2 16059.34 16140.00 16135.82 
Tbed1 1300_3 15606.10 15687.38 15666.23 
Tbed1 1300_4 18143.89 18312.60 18255.10 
Tbed1 1300_5 18191.11 18315.44 18291.05 
Tbed1 i1500 1 26566.69 26824.08 26731.63 
Tbed1 1500_3 28067.68 28362.79 28284.36 
Tbed1 1900 4 28268.36 28518.40 28489.69 
Tbed1 1700_1 37054.60 37751.08 37343.34 
Tbed1 1700_2 34488.56 35076.83 34817.44 
Tbed1 i700_3 34485.24 34977.47 34759.07 
Tbed1 1700_4 38260.98 38860.34 38534.35 
Tbed1 1700_7 21297.30 21437.82 21433.56 
Tbed1 i700_8 20659.96 20823.75 20820.70 
Tbed1 i700_10 22055.41 22274.57 22210.50 
Tbed1 1700_11 17120.15 17189.64 17188.47 
Tbed1 1700_12 18135.97 18232.53 18201.06 
Tbed1 i1000_1 49681.02 50734.33 50104.98 
Tbed1 i1000 2 50842.16 51677.00 51277.80 
Tbed1 i1000 3 47362.62 48141.82 47131.13 
Tbed1 i1000_4 49029. 12 49910.85 49408.86 
Tbed1 i1000 5 50971.44 51824.38 51166.18 
Tbed1 11000_6 27823.84 28051.58 28043.73 
Tbed1 11000_7 27252.32 27521.50 27412.21 
Tbed1 1i1000 8 27373.37 27638.39 27543.01 
Tbed1 i1000 9 26857.09 27127.70 26992.81 
Tbed1 i1000_10 27186.99 27469.49 27242.63 
Tbed1 i1000_11 22180.33 22297.32 22247.61 
Tbed1 11000_12 22160.39 22231.34 22231.18 
Tbed1 11000_13 22657.09 22768.69 22745.61 
Tbed1 i1000_15 22629.44 22706.59 22704.67 
Tbed1 11000_16 21331.81 21389.82 21365.65 


Tbed1 11000_20 21560.86 21618.06 21601.55 


2.4 FARM APT 


作者 进一步 验证 了 算法 机 制 或 参数 设置 对 计算 性 能 的 影响 。 首 先 ， 初 始 解 对 计算 结 

影响 显著 。 作 者 设计 三 种 初始 解 生成 方法 : LH、 线 性 松弛 启发 和 ADD 方法 。 线性 松弛 
启发 方法 是 将 SSCFLP 的 指派 变量 修改 为 连续 变量 ， 求 解 模型 ， 再 修复 多 源 指派 为 单 源 
指派 。ADD 方法 过 程 如 下 : 随机 选择 一 个 设施 ， 进 行 贪心 指派 ， 以 贪心 的 方式 选择 设施 
并 进行 指派 ， 直 到 所 有 需求 得 到 满足 。 比 较 发 现 : 使 用 三 种 方法 整体 差异 不 大 ， 线 性 松 
弛 启 略 优 于 LH 算法 ，LH 算法 优 于 ADD 方法 。 其 次 ， 在 大 邻 域 选择 时 ， 测 试 是 否 使 用 设 
施 排除 机 制 。 测 试 结 果 发 现 ， 两 者 在 求解 质量 和 计算 时 间 方 面 均 无 明显 差异 。 第 三 ， 使 
用 群 解 搜索 取代 单 解 搜 索 ， 即 维持 一 组 精英 解 ， 算 法 随机 选择 一 个 精英 解 作 为 当前 解 进 
行 搜索 。 测 斌 发现 ， 群 解 搜索 略 好 于 单 解 搜索 ， 但 计算 时 间 有 明显 的 增加 。 


3 结论 及 讨论 


本 文 设 计 了 数学 启发 算法 求解 SSCFLP 案例 ,并 使 用 5 个 案例 集 测试 了 算法 性 能 。 算 
法 原理 如 下 : 从 当前 解 中 选择 一 个 局 部 区 域 ， 建 立 该 局 部 区 域 的 子 问 题 模型 ， 使 用 MIP 
优化 器 求解 ， 然 后 更 新 当前 解 。 和 迭代 执行 这 一 过 程 ， 直 到 当前 解难 以 更 新 。272 个 案例 
测试 表明 : 本 文 数学 启发 算法 共 发 现 191 个 案例 的 最 优 解 ， 更 新 36 个 案例 的 已 知 最 好 
解 ; 在 困难 案例 (案例 集 Yang 和 TB4) 和 大 规模 案例 (案例 集 Tbed1) 上 , 解 的 质量 整体 上 最 
优 ， 计 算 时 间 明 显 降低 。 

从 算法 原理 方法 的 角度 ， 本 文 数学 启发 算法 属于 超大 邻 域 搜 索 (VLNS) 算 法 。 本 文 算 
法 与 HMEH 算法 (Tran etal.2017) 类 似 ， 均 是 根据 供需 空间 分 布 确定 一 个 超大 邻 域 ， 在 大 
邻 域 空 间 中 进行 邻 域 搜索 ， 但 也 存在 两 点 不 同 之 处 。 首 先 ， 邻 域 空间 定义 不 同 。HMEH 
算法 选择 的 邻 域 空 间 涉及 15%~20% 的 候选 设施 数量 , 不 考虑 当前 使 用 设施 数量 ; 本 文 算 
法 从 当前 解 中 选择 7 个 左右 的 设施 ， 将 这 些 设施 所 履 盖 的 范围 作为 邻 域 空间 。 对 于 超大 
规模 案例 T4 和 T5，HMEH 算法 定义 的 邻 域 空间 巨大 ， 而 本 文 算法 定义 的 邻 域 空间 相对 
较 小 。 第 二 , 两 者 搜索 方法 不 同 。 HMEH 算法 基于 动态 构建 的 超 图 , 在 并 行 处 理 框 架 中 ， 
使 用 “前 向 选择 ”和 “后 向 选择 ”方法 尝试 改进 当前 解 ， 而 本 文 方法 是 构建 邻 域 空间 内 
的 子 问题 数学 模型 ， 使 用 MIP 求解 器 求解 ， 进 而 更 新 当前 解 。 因 本 文 将 邻 域 空间 大 小 定 
义 在 一 个 可 以 接受 的 范围 内 ， 对 应 的 子 问题 可 以 较为 高 效 地 求解 ， 也 使 算法 效率 较 高 。 
算法 测试 表明 ， 本 文 数 学 启发 算法 在 求解 质量 和 计算 效率 方法 均 优 于 HMEH 算法 。 

SSCFLP 算法 具有 几 个 明显 的 特征 。 第 一 ， 拉 格 朗 日 松弛 方法 在 不 同 案例 集 上 的 求解 
质量 差异 较 大 ， 但 因 其 计算 速度 快 ， 常 作为 复杂 算法 的 初始 解 (Ahuja et al. 2004; Tran et 
al. 2017); oi-aeH In SE SENG Fa te Holmberg et al. 1999; Chen & Ting 2008)。 第 二 ， 除 HMEH 
外 ， 高 性 能 SSCFLP 算法 依赖 于 MIP 优化 器 。 割 平面 (Yang et al. 2012)、 改 进 割 平面 
(Gadegaard et al. 2018)、 核 搜索 (Guastaroba et al. 2014)、 廊 道 方法 (Caserta & Voß 2020) 以 
及 本 文 数 学 启发 均 需 要 MIP 优化 器 。 超 大 邻 域 搜索 有 可 能 较 好 地 求解 SSCFLP,， 但 超大 邻 
域 涉 及 多 个 设施 ， 又 需要 兼顾 几 十 甚至 上 百 个 需求 点 的 指派 ， 组 合 数量 庞大 ， 搜 索 计 算 
复杂 度 极 高 ， 而 MIP 优化 器 在 超大 规模 邻 域 搜索 方面 具有 优势 。 第 三 ， 针 对 SSCFLP 的 
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群 算法 报道 较 少 ， 蚁 群 算法 (Chen & Ting 2008) 表现 较 好 ， 但 该 算 


ue 


Ble; FRAG 


案例 )， 数 学 启发 算法 缺乏 优势 ， 虽 能 获得 几乎 所 有 案例 的 最 优 解 ,但 算法 计算 时 间 明 显 


一 < 


去 未 在 大 规模 案例 上 


行 测试 。 群 算法 较 少 的 可 能 原因 : 区 位 选择 与 需求 指派 的 解 空 间 过 于 庞大 ， 两 个 可 行 
解 差异 很 大 ， 交 叉 操 作 难 以 形成 可 行 解 ， 可 行 性 修复 也 很 困难 ， 从 而 造成 演化 过 程 非常 


群 智能 算法 需要 庞大 的 种 群 才 能 较为 有 效 地 探索 解 空 间 。 
应 到 注意 到 , 本 文 算法 仍 有 需要 完善 的 地 方 。 对 于 小 规模 案例 (如 OR-Lib 和 Holmberg 


` 


WK. 原因 在 于 : 这 些 案例 相对 容易 求解 , 改进 分 枝 定 界 (Holmberg et al.1999)、 核 搜索 、 
HMEH、 割 平面 算法 均 表 现 最 好 ， 直 接 使 用 CPLEX 等 优化 器 也 能 高 效 地 求解 这 些 案例 。 


St 


此 ， 对 于 小 规模 案例 ， 数 学 启发 算法 迭代 求解 MIP 模型 显得 元 余 ， 较 大 幅 地 增加 了 求 
解 时 间 。 未 来 进一步 的 研究 包括 : 分 析 SSCFLP 案例 模型 的 求解 难度 , 判断 直接 模型 求解 
还 是 采用 数学 启发 方法 求解 ; 分 析 邻 域 空间 大 小 对 求解 质量 和 计算 效率 的 影响 ,探索 更 
优 的 邻 域 定义 方法 ， 验证 数学 启发 算法 对 于 SSCFLP 变种 问题 的 有 效 性 ， 如 限制 设施 数 


量 、 要 求 设施 服务 区 空间 连续 等 情形 。 
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A matheuristic algorithm for the single-source capacitated facility 


location problem 


Abstract: This article proposes a matheuristic algorithm for the classical single-source facility 
location problem (SSCFLP). The algorithm starts from an initial solution, and then iteratively 
improves the solution by searching the large-scale neighborhoods. The initial solution is 
generated by a Lagrangian heuristic and the large neighborhoods are explored by solving 
sub-problems. The performance of the algorithm was tested on 5 set of benchmark 
instances. Experimentation showed that the instances could be solved effectively and 
efficiently. For the largest set of instances, the proposed matheuristic algorithm performs 
better than existing methods in terms of the solution quality, the computational time, the 
number of optimal solutions and the number of better solutions. 

Keywords: single-source capacitated facility location problem, matheuristic algorithm; Large- 


scale neighborhood; sub-problem. 
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